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■ Abstract 

m 

We study the possibility that a keV-MeV mass hidden photon (HP), i.e. a hidden sector U(l) 

■ gauge boson, accounts for the observed amount of dark matter. We focus on the case where the 
Qh| hp interacts with the standard model sector only through kinetic mixing with the photon. The 

relic abundance is computed including all relevant plasma effects into the photon's self-energy, 
D I which leads to a resonant yield almost independent of the HP mass. The HP can decay into 

i-C I three photons. Moreover, if light enough it can be copiously produced in stars. Including bounds 

from cosmic photon backgrounds and stellar evolution, we find that the hidden photon can only 
give a subdominant contribution to the dark matter. This negative conclusion may be avoided 
^ I if another production mechanism besides kinetic mixing is operative. 

(N 
m 

■ 1 Introduction 
oo . 

, The evidence for the existence of dark matter (DM) is compelhng[Tl[2]. The various measurements, 
using rotation curves of galaxies, lensing data, cluster dynamics, large scale structure and CMB 
data all agree: about 25% of the energy budget of the Universe is in the form of dark matter, 
and a large part of the dark matter is non-baryonic. Although observational evidence for dark 
\ matter is plentiful, we are still in the dark about its identity. The theory of structure formation 
provides indirect evidence about some of its properties. It favors cold dark matter (CDM) that is 
weakly interacting and non-relativistic at late times. One of the best motivated CDM candidate is 
a WIMP, a weakly interacting massive particle [HIl]. The ligthest super symmetric particle is the 
archetypical WIMP example. WIMPs are in thermal equilibrium with the standard model (SM) 
particles in the early universe. With weak scale masses and interactions WIMPs would have fallen 
out of equilibrium at the right time such that their relic density today is in the right ballpark. 

Although CDM provides a consistent picture of structure formation on large scales, there are 
persistent problems on subgalactic scales. Most notably, CDM predicts too many galactic satellites 
[3 m [71 [8] , a galaxy density profile that is too cuspy [9l [ini [III [E] , and too low angular momenta 
of spiral galaxies [13]. All these problems suggest that CDM may be too cold. This has motivated 
warm dark matter (WDM) [14^ I15j. WDM consists of lighter particles, in the 1-10 keV range, that 
are on the borderline between non-relativistic and relativistic at the time of structure formation. 
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The non-zero velocity dispersion suppresses structure below the Mpc scale. As a result structure 
on sub-galactic scales is damped, resolving the tension between DM matter simulations and ob- 
servations. Note in this respect that even lighter particles would still be relativistic at late times; 
such hot dark matter damps structure on scales much larger than galaxy scales, in conflict with 
observations. 

The standard model of particle physics (SM), while describing collider phenomenology very 
successfully, does not offer a suitable DM candidate, and we must look further. Most SM extensions 
include hidden sectors, i.e. sectors which couple only very weakly, typically gravitationally, to the 
SM fields. Given the importance of symmetries in the SM it is not unlikely that these hidden 
sectors contain gauge groups as well, including some abelian U(l) factors. The usual assumption is 
that such hidden photons (HP) are heavy, and decouple from the standard model particles, thereby 
avoiding all observational constraints. But this needs not be the case. For example, if the hidden 
U(l) is broken by non-perturbative effects, the symmetry breaking scale and thus the photon mass 
is exponentially suppressed, and can be naturally light [1^. At low energies the dominant coupling 
of the hidden photon to the SM will most likely be via kinetic mixing with the SM photon. Kinetic 
mixing is allowed by all symmetries, and can be realized via a renormalizable coupling — thus in 
principle unsuppressed by heavy mass scales — of the form L B xFfiuB'^'^, where and B^j^i, are 
the SM and hidden photon field strengths respectively |17j . 

In this paper we study the possibility that a hidden photon with mass in the keV-MeV range, and 
whose dominant coupling to the standard model is via kinetic mixing, is a dark matter candidate. 
Because "warm" is normally used for thermal relics with masses in the 1 — 10 keV, our keV-MeV 
candidates are more likely "lukewarm". Such hidden photons are too heavy to be measured in 
laboratory experiments; the 5th force searches or the laser experiments designed to search for 
axion like particles are only sensitive to much lighter hidden photons [18^ [T9l I20j . Nevertheless, 
there are astrophysical and cosmological constraints that restrict this scenario considerably. First, 
there is the requirement that hidden photons (HP) do not overdose the universe. In calculating 
the HP relic density we focus on the case where HPs are solely produced through their mixing 
with photons, although we comment on other possibilities. Second, such light and weakly coupling 
particles are produced in the core of stars, and can subsequently escape the star unimpeded, 
providing an extremely efficient cooling mechanism. This alters the evolution of stars, which is 
bounded by observations. These constraints are strong. Indeed, in the whole mass range below 100 
keV the mixing parameter for which interesting (measurable) abundances are obtained is ruled out. 
Third, the HP is an example of decaying dark matter, as it can decay into three visible photons. 
If it is to be the dark matter its lifetime should be larger than the age of the universe, which 
bounds the kinetic mixing from above. And fourth, the decay products of the HP contribute to 
the galactic and cosmic gamma-ray background, again constrained by observations. As we will see, 
all constraints put together practically rule out our hidden photon dark matter scenario, unless 
another production mechanism besides kinetic mixing is operative. 

Under the assumption that kinetic mixing is the dominant production mechanism of hidden 
photons in the early universe, the hidden photon cannot be a cold dark matter candidate with 
MeV mass or above. The reason is that for such heavy HPs the decay channel into an electron- 
positron pair is open. Stability on the scale of the lifetime of the universe can only be obtained 
for an extremely small kinetic mixing. However then the relic abundance is way too small to be 
dark matter. This raises the question whether such super-MeV mass HP can leave any imprint 
on cosmology at all. If the HPs decay after the time of nucleosynthesis and/or matter radiation 
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decoupling (CMB), can they leave measurable traces? Also here we find a negative answer. 

While this work was in progress two related papers appeared. Pospelov et al. [21] also discuss 
the possibility of hidden photons as dark matter. Although there is a large overlap, we extend their 
analysis by including resonance effects (which actually dominate production), and by significantly 
improving the stellar bounds. As a result, our conclusions differ substantially from [21]. In par- 
ticular, whereas they find some hidden photon mass and kinetic mixing parameters for which the 
HP can be the dominant source of dark matter, we rule out this possibility. Chen et al. [22] study 
HPs with masses of ~ 100 GeV as cold dark matter. Their analysis differs from ours in that they 
include additional non-renormalizable interactions, which provide an extra source of production. 
This allows for a long lifetime (requiring small kinetic mixing) while at the same time obtaining a 
large relic density. 

This paper is organized as follows. In the next section we introduce the model. For a proper 
treatment of photon - hidden photon interactions plasma effects should be taken into account and 
the relevant formulas are discussed. In section[3]we calculate the relic density of hidden photons as a 
function of its mass and mixing. We explain that for hidden photons in the keV range production is 
dominated by the resonance regime, which occurs for temperatures such that the thermally induced 
photon mass equals the hidden photon mass. For heavier HPs the main production channel is via 
electron-positron coalescence. In section H] we discuss the cosmological and astrophysical bounds 
on the HP parameter space, including the bounds from overclosure, stellar evolution, and from the 
diffuse gamma-ray background. In section[5]we discuss the effects of non-renormalizable interactions 
between the hidden sector and the standard model. If these additional interactions dominate the 
production of hidden photons parameter space opens up, and hidden photons can be the dark 
matter in the universe. We end with some concluding remarks. 



2 The model 

Consider a hidden sector U(l) gauge boson, a hidden photon, which couples to the standard model 
photon solely via gauge kinetic mixing [T7j with the hypercharge boson. The low-energy effective 
Lagrangian is [23l [24] 

L = -\f,,F^^ - \b,,B^^- + + £^/,2^^^M, (1) 

where F^^y and B^^ are the photon {A'^) and hidden photon {B'^) field strengths. The dimensionless 
mixing parameter sin^o can be generated at an arbitrarily high energy scale and does not suffer 
from any kind of mass suppression from the messenger particles communicating between the visible 
and the hidden sector. This makes it an extremely powerful probe of high scale physics. Typical 
predicted values for Xo realistic string compactifications range between 10~^^ and 10~^ \2h\ [26l 
[271 [28]. 

The most prominent implication of the kinetic mixing term together with the non-zero hidden 
photon mass // is that photons are no longer massless propagation modes. Similar to neutrino 
mixing, the propagation and the interaction eigenstates are misaligned. The kinetic mixing term 
can be removed by a change of basis {A^B} — > {A^,S}, where A^^ = cosXo^ ^-^id 5 = B — smXo-^- 
Since A and A^ differ only by an unobservable charge renormalization we will drop the R subscript 
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from now on. In the {A, S} basis the kinetic term is diagonal but kinetic mixing provides an 
off-diagonal term in the mass-squared matrix, 



l-i^ sin xl ^^gin;,^;^^ cosxo 
/i^sinxo cosxo cos^ Xo 



(2) 



As a result one expects vacuum photon-sterile oscillations |29j as in the case of neutrinos. In the 
following we will use the notation 7 ,7' for the flavor states, which are the quanta of the {A, S}- 
fields respectively. The mass eigenstates are denoted by 71^2, with 71 mostly "photon-like" and 72 
mostly "hidden photon-like" . 

The above discussion applies to the system in vacuum. In the early universe matter effects 
should be taken into account, as the photon is submerged in a thermal plasma. We can include 
the influence of the plasma in the photon propagatioro through the photon's self energy, which 
enters into the above formalism as a complex effective mass. The real part of the mass encodes the 
refraction properties of the plasma and its magnitude is set by the plasma frequency ujp. We will 
refer to it as the plasma mass or the photon mass denoted by m^. For the temperatures of interest 
the main contribution comes from Compton scattering on electrons, which gives 

2 _ / = 4vra(ne/me), (T <C m-e) . . 

"'^ - \ ^ (2/3)a7rr2, (T » me) ^'^> 

where rie is the electron number density. The contributions from photon scattering off other charged 
particles are completely analogous. The full formulas are given in Appendix lAl 

The imaginary part of the photon effective mass is given by 

ujB{uj, T) = -uj {T^iuj, T) - r^{LO, T)) = -uj (e'^/^ - l) r^{LO, T) (4) 



where r^^^\uj,T) is the production(absorption) rate of photons with energy tj in a thermal bath 
at temperature T \30\ I23j. The "damping factor" D(u;,T) may be thought of as a rate parameter 
measuring the effectiveness of the collisions to stop the coherent development of the wave func- 
tion [31]; alternatively it can be interpreted as the rate at which photons would regain thermal 
equilibrium [30]. The relevant definitions are presented in Appendix [Al 

The effective mixing angle in a damping dominated mediuiE@ is given by [23] 

X^(^^^)^xg (^2_,42 + (^D)2 ■ (5) 

It depends implicitly on the energy and the temperature though m.y and D. The imaginary contri- 
bution to the photon mass, wD, is typically smaller than the real parlH, m^, so it only plays a role 
near the resonance = fi where it acts as a cut off. 

Since the plasma frequency is a steep function of temperature, the resonant condition my{Tr) = 
H divides the temperature range into a region of very suppressed mixing angles (m^ » fi for 



^ Throughout this paper we do not discuss longitudinal photons. We expect that including their effects will not 

affect our results significantly. 

^For the small values of Xo ^6 consider in this paper, the mixing is always damping dominated with ujT) 3> ^x^fJ-^ 
^Both the real and imaginary parts come from matrix elements squared, but the real part interferes with the 

identity matrix and therefore involves less powers of the coupling constant. 
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T ^ Tr) and a region which can be considered as vacuum (m^ <^ for T <^ T^). For heavy 
hidden photons with fi > MeV, the resonance happens when electrons are stiU relativistic at a 
temperature = fj,y/3/{2TTa) ~ 8/i, where we used ([3]). For much smaller masses the electrons are 
non-relativistic at the time of resonance, and oc 7"3/2g-me/T_ j-j^g ^j^g exponential decrease 
the thermal mass is extremely sensitive to T. Consequently for a broad range of values for 
in the range 1 - 10^ eV, the resonance happens not far from Tj. ~ 0.2 mg. As we will see, for small 
HP masses the production is dominated by the resonance, whereas for larger masses at slightly 
lower temperatures. This means that in both cases all interesting physics happens at relatively 
high temperatures T > 0.2 rUe, corresponding to the period when electrons and positrons have not 
completely annihilated and their chemical potential is still negligible. 

That the resonance is more prominent for low mass HP masses can be easily understood. 
Considering only Compton scattering the damping factor D ranges from Sira^ /{3m1)ne (at T <C rrie) 
[23] to ~ a'^T'^/{Truj)Log{iTuj/ml) (at T > nie). This gives for the ratio 

which sets the enhancement of the effective mixing angle at the resonance. Since a; ~ T, ([6]) 
shows clearly how the resonance is enhanced at low temperatures, corresponding to small HP 
masses, and losses importance (first linearly, thereafter logarithmically) as electrons become more 
and more relativistic. For our range of hidden photon masses ^ > keV, the resonance temperature 
is bounded > 0.2 mg, and the enhancement factor ([6|) is never extremely large, a factor 10^ 
at most. 

The photon-hidden photon oscillation frequency also becomes modified in the medium, it reads 

wosc = ^\/(/"^ - "^7)^ + • (7) 

Since typically Wqsc ^ many oscillations take place before a photon is absorbed or scattered so 
the two mass eigenstates, which travel at different speeds, are not likely to still have an overlap at 
subsequent interactions. We can then simply treat 71^2 as two different final states. This is true 
except near the resonance, where Wqsc/D is precisely the ratio in ^ times the vacuum mixing angle 
Xo, and therefore is very small unless Xo ^ 2acj/(3me). However, even in this last case we can also 
treat 71^2 as two final states (neglecting coherence) as long as the effective mixing is small because 
the photon- like wave/state/component 71 will be damped much faster than 72. 
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3 Production of hidden photons 



In this section we calculate the abundances of thermally produced hidden photons. We can divide 
the production into three stages. At high temperatures for which m^{T) ^ 11 the hidden pho- 
ton is comparatively massless; photons are very close to being both interaction and propagation 
eigenstates and the effective mixing angle ([5]) is strongly suppressed. Consequently the amount of 
hidden photons produced is negligible small. We can therefore assume that the initial abundance 
of hidden photons is negligible small. This sets the initial condition. As the temperature lowers, so 
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does the plasma frequency, and the system hits the resonance when m^{T) = /x. The photon bath 
present will partly convert into a hidden photon bath, the efficiency of this conversion depending on 
the effective mixing parameter. At lower temperatures, the effective mixing relaxes to the vacuum 
value X ~ Xo) but still the production can be effective, favored by the fact that the expansion of 
the universe is increasingly slower. 

Since we have argued that coherent effects do not play a role, the evolution equation for the 72 
yield is the usual Boltzmann equation one expects from incoherent production, 

dY2 T2 dins ^ 



d\uT H dlnTS 

where we defined Yi^2 = ^1^2/^ the ratio of the hidden photon- like or SM photon-like density to the 
entropy density, and H is the Hubble constant. The dins term on the right hand side incorporates 
the change in the effective plasma degrees of freedom as species decouple from the thermal bath 
[32j . For the broad HP mass range keV-GeV, there is only muon, electron and neutrino decoupling. 
The relevant definitions can be found in Appendix [Bl Note finally that we assumed that 12^1 
always, making the absorption processes of 72s negligible. 

The HP production rate is given by the sum of different contributions 

where n^+,n^- are the number densities of electrons and positrons, cr^jei C72 ! 1^7271 the spin- 
averaged cross sections of the reactions 716 726 (Compton-like production), e+e~ 72 (pair 
coalescence) and e^e~ 7172 (pair annihilation). The brackets denote the proper thermal average, 
which for a reaction a + 6 — > 72 + c is given by [32] 

naUbia^^cV^.^) = j ^^fa j ^^/6 0-72c(s)^^M6l = j dnadnfo (7^2c(s) Wm61 (10) 

with fa^b the Fermi/Bose distribution function = {^e^^-'^/'^ it 1)"-*^ for initial state fermion/boson, 
^^Moi ~ \/ba — — \va X the Mocllcr velocity, Ea^b^Pafi^^afi the energy, three momenta and 
three velocity of the incoming particles, s the center-of-mass energy (not to be confused with 
the entropy density) and gx the internal degrees of freedom of the particle type x. The relevant 
cross sections can be found in Appendix O At temperatures T <C rUe only the scattering process 
contributes to HP production because of the exponentially decreasing electron/positron density. 
At large temperatures however we also have to include the coalescence e^e~ — > 72 and annihilation 
e~^e~ — > 7i72 processes. The coalescence process is much stronger since it is suppressed by a 
smaller power of the fine structure constant a ~ 1/137 but it is only possible for HP masses greater 
than twice of the electron, i.e. for /i > 2me- Therefore in principle we should consider all three 
processed. 

To make the numerics less demanding we can use Maxwell-Boltzmann statistics for the initial 
particles. Leaving aside the resonant production, for which we will develop a simple and general 
formula, this is justified because the production is dominated by the Wien tail of the distribution. 



* We neglect stimulation/blocking factors for final state bosons/fermions since in practice they do not play a 
significant role. 

''The generalization to include other charged particles like muons or pions is straightforward, but for the HP mass 
range of interests can be neglected. 
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The reason is that the production rate is proportional to some inverse power of the temperature 
times a Boltzmann factor e~'^/^, and thus production is dominated by temperatures T < fi. Only 
the particles in the Wien tail of the distributions will have enough energy to produce a HP, and 
their distribution is weh encoded by Maxwell statistics. Note also that in this regime jj, ^ and 
this justifies why we can set = in the cross sections of in Appendix [Cl Under these conditions 
the thermal average ()10p simplifies enormously 

nam. \ 9a9b 1 X, r r f m^ 

^^^"--^ = 1732^W^ Xo " ^^^^^^^^'-^^^^ [-T ) ' ^''^ 

where so = {ma + mft)^ and Ki is the modified Bessel function of the second kind. 

We have computed numerically this average and plotted our results for three representative 
cases in Fig. [TJ For illustrative purposes we have extended our calculations to large values of the 
temperature to cover the resonance as well. Looking at the graph of dy2/dlog^QT we see that 
the resonant contribution grows with respect to the low temperature incoherent part for small HP 
mass. Indeed, where for fi > 2me incoherent scattering dominates production, in the opposite limit 
H < 2me resonant production dominates. For heavy HPs with fj, > 2me the dominant contribution 
to the final HP abundance comes from incoherent production via pair coalescence; the resonance 
contributes only a small fraction. We discuss this case in the next subsection. Below /j, < 2me the 
coalescence process is not possible, and the incoherent production decreases substantially compared 
to the resonance, which is already very peaked and therefore dominates HP production. More 
details are given in subsection 13.21 



3.1 The coalescence region (/i > 2me) 

For heavy HPs with mass fi > 2me pair coalescence dominates. Compton-like and annihilation 
contributions are subleading because of the additional power of a in the cross sections. Moreover, 
even though they contribute to the resonance (unlike coalescence) this contribution remains less 
than ~ 10 %. 

The coalescence cross section contains a Dirac delta of the center-of-mass energy which makes 
the thermal average trivial. We find {n^+ = n^,- = Ug) 

J^K^mJ = ^^K, (I) {^^' + 2m^)V/.2-4m2 . (12) 

Using this expression the the final abundance can be approximated by 



GeV 



A* 



1 dins 



Y2 ~ 1.2 X lO^^x,^ 

\ u / ^/OoffAioff am J \ I \i 

T=T, 



/g^hcs dlnT3 




(13) 



Here we have neglected the temperature dependence of the number of radiation degrees of freedom 
^cff, ^cff (see appendix [B]) and of the electron mass mg. This is because the integral is dominated 
by a narrow interval around T^; ~ /i/3 where the former cannot change dramatically (thus we can 
evaluate them at T^) and the role of rUe in (112p is minimal. To understand this last argument note 
that, due to the thermal component of the electron mass, there is a maximum temperature Tc that 



7 




Figure 1: Differential production rate of liidden photons for three different masses. For /i = lOmg 
(right) the low temperature production through electron-positron coalescence dominates the final 
abundance, the resonance not being extremely peaked. Coalescence happens only for T < 3fi 
because for higher temperatures the electron thermal mass is too high. Production of HPs of 
higher masses proceeds very much in the same way. For ^ = rrie (center) coalescence is never 
possible, and the low temperature production suffers a 0{a) decrease since it is now dominated by 
Compton scattering and pair annihilation. The resonance is more peaked, dominating production. 
Finally for /i = me/3 (left) the resonance has grown very large with respect to the low temperature 
production; the production rate is suppressed by the lack of ambient electrons. 



allows coalescence. Writting m^iT) = m^Q + naT'^l with / C(l) we find 



(14) 



which approaches the limit Tc ~ 3.3/i when /i ^ me^Q. Above this temperature no coalescence is 
possible (this feature is visible in Fig. [1]). Only in the limit T —> Tc is there a significant phase 
space suppression making the square root in (I12p small. Since Tc is much higher than T^, the 
temperature that dominates production, and since the thermal electron mass decreases as T^, we 
do not expect any significant deviation from the RHS of (113^ (unless, of course, /i ~ 2me^o itself). 

As we mentioned before, resonant production happens at > y^l\foL which is always higher 
than Tc- Therefore coalescence is absent during the resonance. Since at the resonance /i = 
this is of course just another way of saying that the photon thermal mass never exceeds twice 
the temperature-dependent electron mass in a thermal plasma, and photon decay is kinematically 
forbidden [331. 



3.2 The resonant region (// < 2me) 



In the small mass regime with /i < 2me electron positron coalescence is kinematically forbidden, 
and thus the biggest contribution to the incoherent production abruptly disappears. For such small 
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Figure 2: Plotted is the function j = as a function of the temperature in units of electron 

mass. The resonant yield I2 is inversely proportional to j, see (fT6]l . 



masses the enhancement factor in ([6]) fjP/ujD ^ 1 implying that the resonance is very peaked, and 
dominates with respect to the incoherent production. In this section we develop a simple formula 
for the resonant production, which is also valid (although not relevant) for higher HP masses. 

As the resonance time is short, we can neglect the change of plasma d.o.f. in the Boltzmann 
equation If the temperature integral is dominated by the resonance it is easier to do it before the 
thermal average of the cross section. Introducing the notation a = y^a with the mixing parameter 
given in ([5]), the solution of dH]) is 



a, 6 



Next we expand the photon mass around the resonance temperature as m^{T) = fi +m^'{T — Tr) + 
The derivative m^' can be expressed as j{T) x m^/T. The function j{T) is plotted in Fig. [2l 
Limiting cases are j(T) = 3 for non-relativistic electrons and 2 in the relativistic case, but notice 
also the spike at intermediate values due to the steep decline in electron density which suppresses 
the resonant production. Evaluating all other quantities at their resonance value the integration 
becomes 



XoV dins /■ , , . 
2<(2) dins 



T=Tr JO 



dT 



(m2')2(r-r,)2 + (L^D)2 



C(3) Hj{T)T dlnTS 



(16) 



T=Tr 



To get the second line we used that he temperature integral just gives ~ -k /{uDm^). This quantity 
only depends on the outgoing hidden photon energy uj, which allows to reorganize the integrals 
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over particle three-momenta 



y [ dn^dn,^av^,, [ ^la^^r^ = / ^^^^ 1 ' = g^!^ (17) 

uD J {271^ LoB J (27r)3 u; eW^ _ i C(3) T ^ ^ 



since T^{e'^^'^ — 1) = D, and we have used = ji <^T ancjfl g-^^ = g-y^ = 2 in the thermal integral. 
Plugging it back in then gives the 2nd expression in (jl6|) . 



The cross section (the production rate) has dropped out of the number density of HPs produced. 
This is a curious result, as the final abundance is independent of the details of the processes involved, 
it only depends on the shape of the resonance. 

In the limit that the electrons are relativistic, i.e. for temperatures T > nig, we can obtain a 
simple formula through = = 2a7rT^/3, 



17^2 ^GeV\ ^2 1 din. 

" V /" / ° V^S^cfrdlnTS ^^^^ 



Y2 « 1.3 X lO^^x,. X 



(18) 



where we set j{T) = 2. Up to kinematical factors this expression is the same as the coalescence 
production yield (jlSp . Note however that the degrees of freedom are evaluated here at the resonant 
temperature and in (jlSp at ~ 0.3//. Since the resonance happens at much higher temperatures 
{Tr- > 8n) where geB,heS are larger, its contribution is suppressed with respect to the coalescence 
yield. 

The resulting relic energy density in hidden photons today can be deduced using the con- 
servation of the entropy per coming volume to get 

= 2.82 X 10^ Q^^2 • (19) 



with h ~ 0.71 the Hubble constant today in units of 100 km/s/Mpc. Since the values of I2 (|13p . (jl8p 
are inversely proportional to the HP mass (up to a small dependence of the ^eflf; ^eff parameters on 
the production temperature), the relic densities we obtain are practically independent of /U. For our 
numerical calculations shown in Fig. [3] we used the full expression ()16p , rather than some limiting 
approximation such as (fTSj) . 



4 Bounds 



In this section we discuss the cosmological and astrophysical bounds on hidden photons with mass 
and mixing parameter (/z, x)- We do not make any additional assumptions about the hidden sector, 
such as the possible existence of millicharged particles or a hidden sector Higgs. 

Note that in our computation of the relic abundances we have neglected possible backreactions 
that convert hidden photons into standard model particles. This remains true until 712 — ni, i.e. 
until 27r^/ieff^2/(45C[3]) is 0(1). Up to a small mass dependence, this happens for Xo ^ 10~^ (see 
the thin dashed line in Fig. [Sj). Of course our computations are not precise outside this regime, 
but we can qualitatively interpret that for such large kinetic mixing the studied interactions keep 
hidden photons in thermal equilibrium with the standard bath, at least for a short period of time. 

® Recall that we do not consider longitudinal modes, which have other dispersion relation. 
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4.1 Decay rates 



The phenomenology is very different for heavy hidden photons {fx > MeV) and lighter ones. This 
has shown up already in the production rate and also affects the decay rate. The mass eigenstate 
72 has a small admixture of the active photon state, and thus can decay. For HP masses smaller 
than two times the electron mass, < 2me, the only decay channel is into three photons through 
an electron loopll|. For larger masses the dominant decay channel is into an electron-positron pair. 
The relevanllfl lifetimes r are given by 

"™ = 11664000vr3 -' « ^'^^ 



m 



e 




r-^.- = (1 + ^1, (M>2™.) (21) 

with Q ~ 1/137 the fine-structure constant. As the HP photon mass increases, new decay channels 
into pairs of higher mass charged particles open up. Decays into heavier leptons are analogous 
to pip and the contribution of decays into hadrons are given by Tj^^ = t~_^^_ x with R{s) 

the experimentally measured ratio a{e^e~ hadrons)/o"(e^e^ — > n~^fi~) as a function of the 
invariant center of energy mass [3l]. Unfortunately, as we will see, such heavy HPs cannot be 
firmly constrained. 



If the hidden photon decays before the onset of big bang nucleosynthesis (BBN), about 0.1s 
after the big bang, they leave no trace in cosmology. For lifetimes between BBN and today (to — 
4.3 X 10^'' s) the decay products can affect BBN and/or the CMB blackbody spectrum. For lifetimes 
longer than the age of the universe the relic HPs contribute to dark matter so their energy density 
cannot exceed the fiducial value f^DM^^ — 0.1. In addition, photons produced in HP decays cannot 
exceed the measured cosmological photon backgrounds. 

In Fig. [3] we plotted in (fi, x)-parameter space the regions for which the HP decays before BBN 
(white upper right part), during BBN (dark gray band), during the time the CMB is unprotected 
to distortions (medium gray band) and after decoupling (light gray band). The region where the 
HP has a lifetime longer than the age of the universe, and thus can (if abundantly produced) be 
~stable dark matter, corresponds to the white lower part of the plot. The decrease in decay time 
as the HP mass passes the pair production threshold is huge. Consequently the most interesting 
region, where the HPs can affect cosmology and/or astrophysics, is the low mass region fx < 2me; 
Let us from now on call HPs in this region "light" HPs. They are the main focus of this work as 
opposed to to "heavy" HPs which essentially evade all constraints. 

^The tt" anomaly also provides a 3 photon decay channel but the amplitude is suppressed by Aqcd instead of mf 
and therefore is subdominant. There is also a decay into two neutrinos, since in general the hidden U(l) should mix 
with the hypercharge U(l). However, in this case the mixing gets an additional suppression factor — A/f )^ 

which for the masses considered is huge. 

® Decays are only kinematically allowed for temperatures below the resonance (T < Tr) such that the photon and 
electron thermal masses are smaller than /i. Since as we will see hidden photon production for T > Tr is negligible, 
this threshold is in practice irrelevant. 
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Figure 3: Bounds on hidden photons in the mass-mixing plane. HPs that reproduce the right 
amount of DM He on the hue Q2h'^ = 0.1 and the region above is excluded by overproduction. 
Above the thin dashed line HPs will interact strongly enough to reach thermal equilibrium with 
the standard bath. The regions labeled Sun and HB are excluded by an excessive HP luminosity 
in the sun respectively in horizontal branch stars in globular clusters. In the region labeled IDPB 
the HP decay products exceed the intergalactic diffuse photon background. This bound assumes 
that the HP relic density is created through the kinetic mixing as discussed in this paper. If 
one assumes other production mechanisms that lead to ^2h'^ = 0.1 independently of x the bound 
extends all down to the light yellow region. We find no bounds above > 2me — 1 MeV. Also 
shown are regions where the HP decay could influence different cosmological epochs: pre-BBN 
(r < 1 sec), BBN (1 sec< r < 3 min) and post BBN (3 min< r < lO^^ sec), CMB-unprotected (lO^^ 
sec< r < 10^^ sec), CMB decoupled until now (10^^ sec< r < 4.3 x 10^'' sec ). See the text for 
details. 
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4.2 Cosmological Bounds 



The relic density cannot exceed the measured dark matter fraction in the universe as measured by 
WMAP [351 [36] which constrains JIcdm^^ ^0.1. Hidden photons that saturate this bound and 
are stable lie in the black solid line in Fig. [3l We see that only light hidden photons can satisfy 
this criterion, and therefore the matter fraction plotted only includes the contribution from the 
resonance (fTSj) . which as discussed before dominates the light HP production. The space above 
this line can be excluded because it leads to too much DM abundance. Strictly speaking this bound 
only applies to hidden photons with the life-time longer than the age of the universe. However, HPs 
decaying after CMB decoupling are also excluded, as they leave their trace in the CMB anisotropics. 
The remaining parameter space above the CMB-decoupling line can be easily excluded as well, but 
this will be better argued after considering possible constraints on heavy HPs. 

Heavy hidden photons decay into electron/positron pairs much before CMB decoupling. One 
might wonder if the decay products can alter the successful picture we have of the post BBN 
cosmology. Let us address this question for HPs of increasing lifetime. 

Heavy hidden photons in the right upper white corner of Fig.[3]have lifetimes smaller than ~ 1.5 
seconds, which correspond in standard cosmology to the moment when the weak charged-current 
reactions that keep protons and neutrons in thermal equilibrium freeze out (at T ~ 0.7 MeV). 
As a consequence of this freeze out the neutron density is fixed (except for the very slow neutron 
decay) instead of decreasing exponentially as it would have done if still in thermal equilibrium with 
protons. Later, nearly all neutrons will convert into ^He nuclei, which provides an indirect test of 
this epoch. This is the earliest time in the history of the universe that we can test so far. Therefore, 
heavy HPs that decay well before do not leave a trace in cosmology. As we have seen, production 
of heavy HPs remains active up to — /i/3. Prom this it follows that HPs with masses below 
~ 2.1 MeV are in thermal equilibrium and thus are present as active degrees of freedom during the 
p-n decoupling (parameter space roughly below the thin dashed line in Fig. ([3])) and they could in 
principle affect it. 

For longer lifetimes the HPs decay during or after BBN acting as extra radiation during the time 
of BBN affecting the BBN yields as well. Moreover, their decay products can spoil the successful 
agreement between observations of primordial abundances and theoretical computations (see [37] 
for a excellent recent review) . The energetic primaries , fi^ , vr^ . . . will initiate electromagnetic 
or hadronic cascades. The former can include photons with energy larger than > 2 MeV which 
will then photodissociate the light elements; the latter can have a variety of effects like destroying 
primordial nuclei, inter-converting protons and neutrons or altering the baryon to photon ratio. 

Usually one derives bounds on the amount of energy released per photon (here ijY2s/n^^) since 
the electromagnetic/hadronic cascades form a universal decay spectrum [381 EH HO] independent 
of the energy of the primaries. This turns out to be an excellent approach for very energetic 
primary particles, but fails for the lowest masses we want to consider of order MeV. For decay 
photons and electrons below the threshold for this universal cascade, but above the threshold 
for photodissociation, the resulting spectrum will be different. The authors of [51] considered 
specific bounds for a late injection of 10 MeV photons from a relic decaying particle and found 
them to be less stringent than for universally cascading photons/electrons. They quote the bound 
E,yn < 10~^GeV for decay photons with energy < 20 MeV from a relic particle whose lifetime 
is > 10^ s and is present with a relic density n relative to the thermal photons. The bounds present 
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r[sec] 

Figure 4: Cosmological constraints on a massive particle of mass /i decaying into SM particles with 
a lifetime r and a relic number density ny. Above the thick solid black line the decay products 
spoil successful nucleosynthesis [38] producing hadronic showers if /i > 1 GeV (below r = 10^ s) 
and electromagnetic cascades if /x > 5 — 30 MeV (above r = 10^ s). The thick short-dashed line 
comes from a more recent calculation [30]. Distortions of the CMB spectrum would be noticeable 
above the long-dashed line [38]. Our predictions for hidden photons for different masses fj, > 2me 
are shown as thin lines. 



in the literature depend on the plasma interactions considered and (slightly) on the experimental 
values of the abundances taken. In Fig. |3]we have reproduced the results of |38^ I40| as an illustration 
and we have plotted the HP energy release per photon ^Y2s/n^^ as a function of the Hidden Photon 
decay lifetime r. 

The region above the thick line below r = 10^ s is excluded for relics with masses > 1 GeV 
which can trigger hadronic cascades via their decay products. The line /i = 1 GeV cuts this line 
for an abundance which corresponds to Xo ~ 10~^^, but this mass is on the borderline of the limit 
of validity of this bounds. A more dedicated analysis is needed to potentially rule out a (small) 
region around these values. Something very similar happens for lifetimes above 10^ s. Here the 
decay products are limited by their capability of initiating electromagnetic cascades which contain 
photons of energy above 2.2 MeV that can destroy deuterium, or above 20 MeV that can destroy 
"^He giving as a product too much Deuterium and '^He. We have also plotted the more recent 
results of [30] as a dashed line. The authors of [38] claim the bounds to be valid for primary 
photons/electrons above 5 — 30 MeV. Again we have a small region which can be be potentially 
excluded, but it is once again at the limit of validity of the bound, and a more carefull analysis is 
needed. 

To summarize, in the high mass region > 2me we find no significant bounds from BBN. 
Basically, if the mixing is large so that appreciable abundances of HPs are produced, the decay 



14 



width is also large and decay occurs before it can cause troubles to our standard picture of BBN. 

Let us now comment on possible distortions of the CMB blackbody spectrum. A perfect black- 
body has fixed energy and number densities of particles which are governed by a single parameter, 
the temperature. The decay products of hidden photons can change the black body spectrum. 
The photon interactions with the ambient electrons and nucleons will tend to bring them back to 
their blackbody values, although at a different temperature. Compton scattering is effective in 
redistributing energy until the universe is ~ 10^ seconds old, but cannot change the photon num- 
ber. The interactions responsible for this are bremsstrahlung e~p^ — > e~p^j and double Compton 
scattering e"7 £"77, and they become ineffective at around 10^ seconds at ~ keV temperatures. 
Therefore, if the HP decays before 10^ seconds the blackbody will be eventually reestablished. If it 
decays between 10^ and 10^ seconds the photon spectrum will regain a thermal shape, but with a 
non-zero chemical potential. Finally if the decay is later than 10^ seconds there will be a noticeable 
non-thermal distortion in the CMB called "Comptonization" or y-distortion [52]. These last two 
possibilities can be strongly constrained by the precise determinations of the blackbody spectrum 
of the CMB by the FIRAS spectrometer on board of the COBE satellite [53]. The bounds on the 
energy injected were derived in [38] and are also shown in Fig. [5] as a long dashed line. We again 
see that heavy hidden photons cannot produce noticeable distortions. 

Let us now examine BBN and CMB bounds for light HPs (/i < 2me). They have masses 
below the electron mass, so they cannot produce photons with energy greater than 2 MeV which 
could photodissociate deuterium. In Fig. [3] we already see that light HPs that decay before the 
CMB decoupling lie above the thermalization line, so during their resonant production at > 8/U 
they reach thermal equilibrium with photons. This fixes their abundance rather than (jlSp . Later, 
when they decay, their energy density is at least as large as that of the photons. Using ^ ~ mg 
and and examining the CMB bounds of Fig. |4] it becomes clear than the region labeled 

FIRAS in Fig. [3] is completely excluded. Light hidden photons with larger kinetic mixing will decay 
when photon-number interactions are active and therefore will not distort the spectrum. However, 
if they decay after BBN (above the white dashed line) they will certainly distort standard BBN 
acting as invisible degrees of freedom. Since they are in thermal contact until T ~ 8/i ~ 5 MeV 
and then suddenly decouple they miss the injection of entropy of electron positron annihilation, 
and they contribute to the energy density during BBN as 8/7 that of a neutrino species. At least, 
this would be the case if they were massless. Kolb and Scherrer showed [H] that a neutrino with 
mass in the range 0.1 — 10 MeV affects BBN more strongly than a massless neutrino, and the same 
conclusion remains valid for our light HPs. Similar conclusions can probably be drawn even if 
they decay before BBN although a detailed study would be required for this bounds to be robust. 
Finally if Xo ~ 10~^ ~ 10~^ they would decay immediately after their resonant production, before 
p-n decoupling leaving no trace in BBN. 

4.3 Gamma Ray Bounds 

The decay of light dark matter hidden photons 72 — > 37 would contribute to the observed gamma- 
ray background. Bounds on decaying dark matter have been extensively discussed in the litera- 
ture [451 [56l [571 [58l [59l [50l I51j . Most of these works concentrate on 2-body decay. Line searches 
in the spectrum are done to bound the decay rate and mass of the decaying DM particle. Hidden 
photons decay into three photons, and the resulting spectrum has no sharp line feature. Instead of 
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applying the line searches to our 3-body decay spectrum, we demand that the gamma ray flux from 
decaying HPs does not exceed the total observed gamma-ray flu^H. Model-independent bounds 
were calculated in [52] imposing that the photon flux from the DM decays does not exceed the 
measured intergalactic diffuse photon background (IDPB): 



for a particle of mass ^, lifetime r and lo the energy of the decay photon. The energy spectrum of 
the decay photons is peaked around lo ~ ^/3. Ref. [52] assumed Jl/i^ = 0.1, setting the last factor 
in (j22p to unity. The used energy bin width is large enough so that the bounds for 2-body and 
3-body decay are nearly the same. 

In Fig. [3] we plotted the curve in (//, x)-plane for which the HP decay is 10% and 100% of the 
total observed flux. The parameter space above this line is excluded. In our calculation we did 
not assume that HPs are all of the dark matter setting $72^^ = 0.1 in the RHS of (I22p . but instead 
calculated the relic abundance as a function of (fJ-,x) from our calculations of section\^ 

4.4 Astrophysical Limits 

The production of weakly interacting particles in stellar interiors can substantially affect stellar 
evolution [53^ [5l] . Hidden photons with small x will be scarcely produced in the dense plasmas of 
the stellar interiors, but they will leave the star unimpeded contributing directly to the star's overall 
luminosity. On the other hand, only photons of the photosphere (and neutrinos) can contribute to 
the standard energy loss. Therefore, naively, the hidden photon luminosity is enhanced at least by 
a volume/surface factor and a further (pmside//0surface)"(7inside/rsurfacc)"* (with p a typical particle 
density and n, m > 1) with respect to the photon luminosity. 

The effect of such an energy loss in the case of the Sun was studied in [23j . In a main sequence 
star like our Sun the existence of an exotic luminosity simply means that the fusion of hydrogen 
into helium is proceeding faster than required to account only for the observed photon luminosity. 
It is well known that the Sun has been shining for more than 4 billion years, and this constraints 
the rate of helium production (since there is still some H in the Sun). This translates into a bound 
Cry < Cry for the hidden photon luminosity. Integrating over a solar model [55] gives the bound 
labeled Sun in Fig. [3l An improvement of roughly one order of magnitude in the bound can be 
achieved by matching the photon luminosity with the neutrino flux, also proportional to the rate 
of hydrogen fusion, as done for the axion case in [56]. However we already see that the bound 
degrades very fast for ^ Tq ~ 1.3 keV, leaving space for hidden photon dark matter. To close 
the window we need to consider hotter stellar interiors. 

We shall focus on horizontal branch (HB) stars in globular clusters which provide the strongest 
limits for energy loss for intermediate temperatures. The core of a HB star burns helium into 
carbon and oxygen at a quite constant temperature of 8.6 keV and density ~ 10^ gr cm~^. With 
these characteristics a HB core is still a classical plasma and the Compton process provides the 

^The bound obtained is a factor 10~^ weaker than that obtained for line searches for 2-body decay. A similar 
suppression factor is obtained when we try to adopt the line-search result to 3-body decay, as the maximum hits per 
bin is reduced by a factor lO"'^ with respect to the 2-body decay. 




(22) 
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most important HP production mechanism. At such smah temperatures the thermahy averaged 
production cross section can be computed precisely taking into account the correct statistics for the 
photon. The invariant mass is practicahy independent of the electron momentum s ~ + 2m(,oj 
so J drig factorizes out. Also the relative velocity is simply v^^^ ~ 1. We find for the energy loss 
per unit mass 

d^ = ^^i (23) 

which has to be smaller than 10 erg g^"*^ s^"*^ [53j. Note that in a HB core the plasma frequency 
at the core is (j ~ 2 keV so we can take Xcfr — X for masses /x ^ 2 keV. The resulting bound is 
plotted in Fig. El 

Usually hotter stars are also more dense and, having larger plasma frequencies, it is possible to 
have a resonant production for ^ > 2 keV. This can be the case for red giants before He ignition 
(T ~ 8.6 keV, a; ~ 18 keV) and white dwarfs (T ~ 1 keV, u; ~ 23 keV). All these stars have 
cores supported by electron (or neutron) degeneracy pressure. Photon interaction rates in these 
environments are suppressed by Pauli blocking factors of fermions, making the HP energy loss 
calculation more involved. Moreover, the HB bound constraints HPs as a sizable fraction of dark 
matter for ji < 100 keV. Only in supernovae cores (up to T ~ 10 MeV, a; ~ 1 MeV) can the density 
be so high to get resonant production of such massive HPs. Unfortunately, in this case we find that 
the emission of HPs cannot compete with neutrinos, and the bounds are not competitive with the 
IDPB bound described above [21j . 



5 Non-renormalizable operators 



In our treatment so far we assumed that kinetic mixing is the only relevant coupling between the 
standard model and the hidden sector. Kinetic mixing can be generated at arbitrary high energy 
scale by "messenger fields" which are charged under both the SM U(l) and hidden sector U(l). 
However, these same messenger fields generate at one-loop order also a non-renormalizable term 
in the Lagrangian. It is not impossible that these additional interactions dominate the production 
of hidden photons. If so, this opens up parameter space and allows for the possibility of HPs with 
mixing parameter smaller than x ~ lO"'^^ — 10~^^ as warm dark matter (hidden photons with 
larger mixing are still excluded as they overdose the Universe). Let's discuss this scenario in a bit 
more detail. 

For definiteness, consider two messenger fields with masses M, M' and charges (1, 1) and (1, —1) 
under weak hypercharge and the hidden U(l) respectively. The kinetic mixing generated at one 
loop is then [T7j 

In the second step we assumed degenerate masses 5M^ = W'^-M'^ <C M^. In the last step we took 
the hidden sector gauge coupling g' to be of the same order as the hypercharge coupling gy, with 
1/60 = a{Mz) < a < q;(GUT) < 1/25. To get small enough kinetic mixing x ^ 10"^^ requires 
very degenerate messenger fields 5M/M < 10"^, very small g' , or a combination of both. 

The same messenger fields also give rise to non-renormalizable operators in the Lagrangian, for 
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example (up to order one factors) 



A significant amount of hidden photons can be produced through the above interaction provided the 
reheat temperature of the universe is high enough. Since the HP production rate ~ aa'{T/M)^T 
depends on T to the ninth power, for temperatures T <^ M these processes are neghgible, and ki- 
netic mixing dominates the production of HPs. The non-renormaUzable interactions are in thermal 
equilibrium in the early universe if the interaction rate exceeds the Hubble rate H ~ T^/rrip. This 
happens for 

that is for reheat temperatures Trh ^ O.IM. If equilibrium is established the number density of 
hidden and SM photons is of the same order of magnitude at the time the reactions (f25]l freeze out. 
Rewriting to entropy density, and using (|19|) this translates into a HP relic density of 

with hcs ^ 200 the effective entropy degrees of freedom (see Appendix [B]) at freeze-out. The dark 
matter density Q2h'^ = 0.1 is obtained for masses in the ^ ~ keV range. 

For larger HP masses, the right relic density can only be obtained if production occurs out of 
equilibrium. Production is most effective at the largest temperature which is the reheat tempera- 
ture. The estimated relic density is [22] 



1-3/4 / 



(fi.n»„o,~6xlO-VHifvi(^)'(lf) <^«) 

which can easily be of order 0.1 for a wide range of masses, depending on the parameters g' , TrH) M. 
Without an UV completion of the theory no definite predictions can be made. The result is 
independent of the mixing parameter. From (j24|) it can be seen that the mixing parameter depends 
on the unknown mass degeneracy 5M/M of the messenger fields. 

The constraints on this scenario are the same as when kinetic mixing is the only production 
source. As mentioned before, the mixing parameter should be smaller than x ^ 10~^^ — 10"^^ to 
avoid overclosure. The resulting HPs have a lifetime larger than the universe, and can be the dark 
matter. The stellar constraints from the sun and horizontal branch stars are unchanged, but the 
gamma ray bounds are slightly different than before. In Fig. [3] we used the HP parameters 
to calculate the relic density produced via kinetic mixing, which in turn was used to determine the 
gamma ray bound for the parameters at hand. In the present scenario the relic density gets an 
additional and unknown contribution from early production via the non-renormalizable interactions. 
Therefore in Fig. [3] we plotted the gamma ray bounds taking a fixed value of the relic density 
^2h^ = 0.1, as well as the bounds when Q.2h^ is computed from production via kinetic mixing and 
thus depends on {XitA- Parameters above the gamma-ray bound for which ^2/1^ is dominated by 
kinetic mixing are still excluded. Parameters in between the two lines are excluded as the dark 
matter candidate, but for smaller relic densities the gamma ray bounds are still evaded. Finally, 
the parameters below both bounds are unconstrained by the 7-background. 
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The discussion so far concerned hidden photons with masses /z < 2me which could be the dark 
matter. Of course, the production of heavier HP can also be enhanced by the non-renormalizable 
interactions. The resulting constraints from the CMB and BBN are very model dependent, as they 
depend on the unknown parameters M, 5M, Trh- The constraints for a fixed number density can 
be found in the literature [3811391 HOl [37] . 



6 Conclusions 

In this paper we studied the production and cosmology for hidden photons with a mass in the keV- 
MeV range. The main motivation for the undertaking of this project is that such hidden photons 
can be the lukewarm dark matter of the universe. 

In our study we focused on the case that the main interaction between the hidden sector and the 
SM is via kinetic mixing with the photons. In particular the relic abundance is computed under the 
assumption that kinetic mixing is the dominant production mechanism. For a proper calculation 
all relevant plasma effects to the photon's self-energy have to be taken into account. The net result 
is that for sub-MeV hidden photons the production is dominated by the resonance which occurs at 
temperatures such that the plasma mass of the photon equals the hidden photon mass. The final 
abundance is almost independent of the HP mass. The dark matter abundance can be obtained 
for mixing angles of order x ~ 10~^^ — 10~^^. For heavier, super-MeV hidden photons, the main 
production channel is via electron-positron coalescence. 

Hidden photons are an example of decaying dark matter. Light HPs can decay into 3 photons 
via the kinetic mixing interaction; heavier HPs with a mass /i > 2?ne can also decay into a electron- 
positron pair with a relatively large decay rate. For the HPs to be the dark matter its lifetime 
should be of the order of the age of the universe or larger. This excludes the heavier HPs with super- 
MeV masses. What is even more, we showed that such heavy HPs leave practically no imprint on 
cosmology even if they decay after the time of nucleosynthesis and/or decoupling. This is different 
for sub-Mev hidden photons. Indeed, light HPs with mixing parameter x ~ 10~^^ — 10"^^ are 
stable enough to be dark matter. But they do decay, and the decay photons contribute to the 
diffuse gamma-ray background. As this background is measured in the relevant energy range, this 
puts bounds on the HP relic density, and thus on the mass and mixing parameter. We find that 
HP with masses ^ > 100 keV are excluded as dark matter. 

Apart from the bounds on the relic density, the lifetime of the HPs, and the diffuse gamma ray 
background there are in addition bounds from stellar evolution. If light enough, HPs are copiously 
produced inside stars. As their interaction rate is small, they subsequently leave the star (almost) 
unimpeded. Thus HPs present an effective cooling mechanism for stars, which is bounded by 
observations. For our study in particular the bounds from the Sun and horizontal branch stars in 
globular clusters are important; they exclude all parameter space for masses ^ < 100 keV. 

All bounds put together exclude HPs as main source of dark matter. There is a loophole in this 
argument though. We assumed that production is dominated by kinetic mixing. But if the reheat 
temperature of the universe is high enough non-renormalizable operators, which are always present, 
can play a role too. In particular, if HPs are pre-dominantly produced via these non-renormalizable 
interactions, parameter space opens up, and HP can still be the dark matter in the universe. Even 
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if this is not the case, HPs can provide a fraction of the dark matter, which can be testable in the 
future through its 3-photon decay. 
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A Photon self energy in the primordial plasma 

The photon self-energy in the early universe plasma receives its lowest order contribution in a from 
Compton scattering off the thermal electrons |571I58| . It turns out to be different for transverse (ttt) 
and longitudinal excitations (ttl). In this paper we only focus on the first case. The computation of 
the real part is nicely reviewed in [53]. One finds that it is closely related to the plasma frequency 

Re TTT = ~ i^l ^1 + {ylk^lup-)^ (29) 
cvl = — r (/e+ + fe- ) (v - Iv') ede (30) 



with k, uj the photon momentum and energy, e the electron energy, /g± the Fermi-Dirac distributions 
of electrons and positrons and the velocity is = \/l — m^/e^. Here v^, is a typical electron 
velocity given by = LOi/up and tui and G{x) are defined through 



col 



^ rife++fe-)(lv'-v')ede (31) 

(32) 



2x 1 — X /! + y/x 
1 - ^ - TT^Log 



G(x) = - . ^ , 

X [ 3 2Vx ^yi- y/x^ 

Note that G{x) ranges from to 1, so in practice the effective mass is never far from the plasma 
frequency. In the cases we are interested in photons have energies much larger than cop and therefore 
one can take A; ~ lij in (I29p . The effective mass is then effectively independent of oj. 

The imaginary part is given by (HI) though the difference of the photon production and 
absorption rates F^. Considering the reactions {i} — > {/} + 7 where {i},{j} are sets of possible 
initial/final states they are given by 

r^iu;,T) = ^Y.Il I d%l-A^(W-{/} + 7)P I n f n i^^fb)] (33) 
{i} {/} / \be{f} / 

r^('^,r)=^5^^ / dn^fiMi + if} ^ I n I n(i±^)) (34) 

{/} {i} Vae{/} / \be{i} J 

(35) 
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where A4 are spin-averaged matrix elements, fx are either Fermi/Bose distributions, the it is to be 
chosen to account for the stimulated boson emission (+) or for the fermion blocked emission (— ). 
The phase space integration is given by 

a€{i} 6e{/} \ae{^} ^ ^ "J \be{f} ^ ^ "J 

where Px{Px) and are the 4(3) momenta and energy of the partice species x and gx the number 
of internal degrees of freedom. 



B Degrees of freedom in the early universe 



The expansion rate is governed by the Hubble parameter, which neglecting curvature effects is 
given by 

= (37) 

where is Newton's constant and p is the universe energy density. Usually one expresses the 
total energy and entropy densities p, s by normalizing them with the contribution of a thermalized 
relativistic bosonic species getting the so-called effective number of species 

5«ff(^) = ^ P ; '^efT(r) = s . (38) 

The factor 

dLogs _ I T dfeeff , . 

dLogTS 3/ieff dT ^ ^ 

is different from 1 only near the annihilation thresholds T ~ nix of standard model particles. 



C Cross sections for hidden photon production 



Hidden photons interact though their kinetic mixing with photons with all standard model charged 
particles, i.e. though the electromagnetic current j^™ as 

Therefore for every SM process in which a photon is emitted, there is also a HPs production channel 
with a coupling reduced by x{'^-,T) given by ([5]). 

The cross section of electron-positron coalescence into a HP is given by 



cj^,(s) = 4^2^ - 4mf [l + ^ j 6{s - p") (41) 
with 5{s — M^) a Dirac delta function. Recall x is a function of uj which in this case is = ^/s. 
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The cross section of e^e annihilation into a (massless) photon and a HP is given by 



a. 



7271 (^) = ^ I 2 ^,..9^ ^Log 



s(s — 4m2) 



2mp 



s(s + 4m2) + 



S\J s(^s — 4ml) 
(42) 

where x = xi^,T) from dS} and uj = {s + fi'^)/{2^). If the HP mass is big enough (/x > 2me) 
then we find t and u-channel divergences at the threshold s = fj?, reflected in the denominator 
of the first factor on the right hand side. In a thermal bath, these divergences are cut off by the 
photon thermal mass. Since the main contribution to the non-resonant HP production comes from 
temperatures smaller than the resonance /i ^ m^ and the inclusion of m-y affects the cross section 
mainly in the divergent electron propagator, which we can correct applying the prescription 



s — n + 2m^iJ, 



(43) 



to the first term. 



Finally, the cross section of Compton-like production of hidden photons is given by the more 
complicated formula 



C7^2e(s) 



2^,2 



2vra^X 



(— [s^ + I5s'^ml - sml + + /i^ (7s^ + 2sm^ - mg)) + 

s(l + (3)+ml- 



(s — m2)3 

+2 (s^ — 6smg — 3mg — 2/i^(s — ml — /i^)) Log 



2meVs 



(44) 
(45) 



where (3 is an invariant defined by s^0^ = {s — [m^ + /i)^)(s — (me — ^)^), which is proportional to 
the HP momentum. In this case the HP energy is ti; = (s + — ml.)/{2^/s). 
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